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ABSTRACT 

In this paper first the plane elasticity problem for an orthotropic 
strip containing a crack parallel to its boundaries is considered. The 
problem is formulated under general mixed mode loading conditions. It 
is shown that the stress intensity factors depend on two dimensionless 
orthotropic constants only. For the crack problem the results are given 
for a single crack and two collinear cracks. The calculated results 
show that of the two orthotropic constants the influence of the stiffness 
ratio on the stress intensity factors is much more significant than that 
of the shear parameter. The problem of loading the strip by a rigid rec- 
tangular wedge is then considered. It is found that for relatively small 
wedge lengths continuous contact is maintained along the wedge-strip 
interface, at a certain critical wedge length the separation starts at 
the midsection of the wedge, and the length of the separation zone 
increases rapidly with increasing wedge length. 


1. Introduction 

From the plane elasticity solution of an infinite orthotropic 
medium containing a series of collinear cracks it is known that the 
material orthotropy has no effect on the stress intensity factors [1] 

(see also [2] and [3] for the case of a single crack). On the other hand, 
the solution given for an infinite strip has shown that if the medium is 
bounded then the stress intensity factors would be influenced by the 
material orthotropy [4,5]. If the (collinear) cracks are perpendicular 
to the boundaries of the strip and if the external loads are independent 
of the coordinate parallel to the boundaries, then the stress intensity 
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factors and indeed the stress component in the plane of the cracks turn 
out to be invariant with regard to a 90-degree material rotation [4]. 
This result is not valid if the strip is loaded along its boundaries 
(e.g., bending of the strip by transverse shear loads) [5]. Even though 
there exists a general solution for an orthotropic strip containing an 
inclined crack [6], a systematic solution of an orthotropic strip con- 
taining cracks parallel to its boundaries investigating such factors as 
the relative size and location of the crack, the types of loading, and 
the material constants does not seem to exist in literature. Such a 
problem would have applications to, for example, the study of delamina- 
tion fracture in layered composites, wood and other homogeneous ortho- 
tropic solids. 

In this paper, first the general problem of an orthotropic strip 
containing a crack located at an arbitrary distance from and oriented 
parallel to the boundaries is considered. The effect of interaction 
of more than one crack on the stress intenstiy. factors is studied by 
considering two collinear cracks in the strip. The important problem 
of splitting the strip by a rigid wedge is then investigated. For a 
wedge of rectangular cross-section which is short relative to the crack 
length it is shown that the wedge remains in contact with the strip 
along its entire length. However, for wedges longer than a critical 
length the contact is shown to be restricted near the ends of the wedge 
only. It is also shown that neither the critical wedge length nor the 
length of the contact region is dependent on the wedge thickness, 

2. Formulation of the Crack Problem 

Defining a stress function ^ by 

_ 3^^ _ 3^(^ _ 3^(^ /^\ 

'^ll ■ 3V * *^22 ■ » *^12 ■ ■ 3 Xj3X2 * ^ ‘ ^ 

the compatibility condition for the generalized plane stress problem 
for an orthotropic plate may be expressed as [7] 
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Egg 3Xi 


+ ( 
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3^ 


3X,^3X. 


3‘♦(^ ^ 
3X2^* 


= 0 


where the notation for the engineering material constants is given by 
the following stress-strain relations: 


^11 " ^'^11 ■, ''12^^22 " '’13°33^ » ••• * 

9-1 

" G,o '^IZ ’ ••• 


C3) 


For the plane stress problem introducing [8] 


E = »^i-iEg2 , V = /v^gVg^ , s'* = Eii/E.22 



and defining 

Xi 

^ ^ ^ ^2*^ * °xx " ‘^ll^^’ % " '^'^22 ’ °xy " ^^12 

from Cl) and (2) it follows that 


(4) 

(5) 



xy 3X3y ’ 


(6) 


1)^ ■*■ 2'" 35^^ ° • (7) 

It may easily be verified that with the transformation (5) equations (6) 
and (7) remain valid for the plane strain case provided the stiffness 
ratio 5 and the shear parameter < are defined as follows [8]: 


54 = ^‘^23^32 

^2 ^"''13''31 * 


K 


_ 1 r ^11^22 j, 

2 ‘•(l-v^3V3T)(l-V23V3g)-‘ 


Va 


/ 1 ^21~^23^31 ^12~*^13^32x 

'G ” F ” F ■ 1 • 

h2 42 4l 

( 8 ) 
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In this case the effective stiffness E and the effective Poisson's ratio 
V are defined by , 


E = [■ 


^11^22 






(9) 


Also introducing the new displacements by 
u = , V = U2/»^ » 


( 10 ) 


The strains and the stress-strain relations may be written as follows; 


'XX 


- ly. 

"11 “ 3x ’ ®yy 


_ ^22 _ 3v 


3y ’ xy 12 ^3y 3x 


m|^ + |J)/2, (11) 


^xx " E ^‘’xx"'’v' ’ ^yy ' E t^yy'^xx'’ ^x 


_ (k+v) 


xy E ’ 


( 12 ) 


For the strip shown in Fig. la in terms of standard Fourier inte- 
grals the solution of the differential equation (7) may be expressed as 


<l>(x,y) = (f^(x,y) + <|.2(x,y) 

00 

<{*l(x,y) = ^ I Cj(a)e“^^j"'*“^ do , -H2<y/vT<Hp 
*00 

-|a|ySj“iax 


♦2(x.y) = i I Aj(«), 


do , y > 0 , 


1 f” 2 lajys.-iox 

<l>2(x5y) = 17 j .2 do , y < 0 

-eo 0"^ 


(13) 

(14) 


(15) 
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where are the roots of 

s'* - 2 ks2 + 1 = 0 (16) 

and are given by 

"^2 

S^ = (k+i/k 2-1 ) , Sg = (k-/k2_1) , S 2 = -S^, S^ = -S2 . (17) 

The functions (j>^ and respectively associated with a strip without 

any cracks and an infinite plane containing collinear cracks along 
the x-axis. The unknown functions Cj and ^y (j=l,...,4) are determined 
from the following eight boundary and continuity conditions (see Fig. la 
and equations (5)) 

*^22^^1’^1^ ~ ~ ®22^^1*"*^2^ ~ 0i2(x^»“H2) “ 0, 

-« < x-j < « , (18) 

<722(^1*'^^) ” 022(^1*”^) » ^12(^1*'^®) “ *Jl2(^l»“^^» ““ ^ ^1 ^ 

(19) 

022^^1*^^ “ o(x-j) t ^12(^1*®) “ T(x-j) » -a < x-j < a , (20) 

U2(^i »'^0) “U2 (xi »“0) “ Oj u-j (x-j ^+0)”U^ (x-j »“0) — 0 , |x-jj > a , (21) 

The six homogeneous conditions (18) and (19) may be used to elim- 
inate six of the eight unknown functions. The mixed boundary conditions 
(20) and (21) would then give a system of dual integral equations to 
determine the remaining two functions. By defining the new unknown 
functions 

G2(xy) = -^ Cu^(Xt.+0)-u^(x^,-0)], G^(x^) = [u2(x^ ,+0)-U2(x^ ,-0)] 
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the problem may also be reduced to a system of integral equations in 
G-j and 62- In this case it is seen that (21) is equivalent to 

r® 

Gj(x^) = 0, Ix-jl > a , J Gj.(x^)dx-j = 0 , (j = 1,2) , C23) 

• -a 

and (20) gives the desired integral equations. 

The displacements and stresses contributed by <j>.j are continuous 
on X2 = 0 plane. Thus, the functions Ag and are eliminated by sub- 
stituting from (15) and (17) into the continuity conditions 09)» giving 

Aj “ C(Si+S2)A-j + 2 s 2A2 ]/(s 2 ~Si ) , 

A4 = [2s^A^ + (s^+S 2)A2 ]/(s^-S 2) . (24) 

Substituting now from (6) and 05) through the Hooke's law into (21) 
and by using (22), A^ and A2 may be obtained in terms of the Fourier 
transforms of G.j and G2. The functions Cp...,C^ are then determined 
from the boundary conditions (18), again in terms of the Fourier trans- 
forms of G-j and G2. Noting that G-j and 62 are zero for lx-j|>a, the 
expressions found for A-j, A2, C-j,...,C^ are of the following form 

2 f® • * 

Aj(a) = Z 2j.^(a) J G^(t)e’“^ dt,j = l,2, 

-a 

a 

Cj(a) = Cj^(a) I G^(t)e^“^ dt, j = 1 4 , (25) 

•3 

where ajj^ and Cjj^ are known functions. Substituting from (5), (6), (13)- 
(15) and (25) into the remaining boundary conditions (19) and performing 
a standard asymptotic analysis to separate the singular parts of the 
kernels we obtain the following systsn of integral equations to determine 
the unknown functions G-j and G2: 
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7 J ^ t,-x, ' * ■*■ 

-a 

■ Emi’ , -a < x^ < a , 

1 f^ 5G„(t, ) 

7 j ^ t^-x^ '*’ ‘‘' 


1 

Em2 



-a < x-j < a , 


(26) 


where k.^ ((i,j) = (1,2)) are known bounded functions and the constants 
m-| and are given in terms of the roots s-j, s^ (see (17)). The deriva- 
tions leading to (26) as well as the expressions for the kernels k.. are 

I J 

rather lengthy but very straightforward. Therefore, they will not be 
reproduced in this paper. The details of the analysis may be found in 
[9]. 

One should emphasize that the kernels k.. are dependent on the 

* V 

shear parameter < through the characteristic roots s-j and S 2 » the stiff- 
ness parameter 6, and the dimensions a, H-j and H 2 , (H-j+H 2 =H). However, 
the analysis shows that aside frtjm k the remaining variables enter into 
the expressions of the kernels through the ratios H 1 /H 2 and 6H/a only. 

The solution of the system of singular integral equations is of the 
following form: 

f^i(x-|) 

G..(x,) = ^—4: , j = 1,2,... (27) 

^ ^ (a2-x^2)* 

where and F 2 are bounded unknown functions. After normalizing the 
interval (-a, a) (26) may be solved numerically by using a Gaussian 
quadrature formula [10]. 
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3. Stress Intensity Factors 

The Modes I and II stress intensity factors at, for example, the 
crack tip X] = a are defined by 

k-jCa) = lim /2(x,-a) a 22 (x-,, 0 ) , 

X-j-^ 

k 2 (a) = lim i/2{x,-a) a-j 2 (x, ,0) . (28) 

x-j->a 


Observing that equations (26) give the stress components a22(x-j,0) and 
cti 2 (xi, 0 ) on the plane of the crack for |x^| > a as well as [x^j < 0, 
substituting from (27) into (26) a simple asymptotic analysis would 
show that 




Em^ F^(-a) 


^12^^!’^^ “ Em25[ 


✓2a(-x-j-a) 

p2(-a) 

/2a(-x^-a) 


F^(a) 

/2a (x^ -a) 

_!2<!L 

*^2a(x^-a) 


+ 0(i/xi^-a^) ] , 
+ 0(/x.j^-a^)]. 


(29) 


Thus, from (28) and (29) we obtain 


Em 


k-j(-a) = F^(-a) , k^(a) = - 

5/a 


Em-| 

5/a 


F,(a) 


k2(-a) 


Em^s Em«6 

—— F2(-a) , k2(a) = ^ F2(a) . 


(30) 


From (22), (27), and (30) the crack opening displacements in the 
neighborhood of the crack tip, for example, x-j = a may be expressed as 
follows: 
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(31) 


U2(xi,+0)-U2(x^,-0) = 
Ui(x^,+0)-u^(xp-^) = 


5kr(a) \ 

— /eu-x^) + 0((a-x^) 

kpCa) 3/2 

^.TO^+0((a^Xi) 


) . 
) . 


If U is the work done by the external forces and V is the internal 
potential, for a crack extension da the crack closure energy may be 
expressed as 


a+da 


d(U-V) = I ^ { 022(^1 


+ a-j 2 (x-j, 0 )[ui(x-j-da,+ 0 ) - u^(x^-da,- )]}dx-| 


(32) 


From (31) and (32) the rate of internally released or externally added 
energy (for the extension of the crack at one tip) may then be expressed 
as follows: 


G = 


_d. 

da 


(U-V) 


6k , 2 k«2 


7E 


m. 


1 


6ni2 


(33) 


In the isotropic materials 5=1, k = 1, m^ = m 2 - 1/4, and for the 
plane stress case under consideration G becomes 


G = f (k^2 + k2^) 


(34) 


4. Collinar Cracks 

In formulating the problem, aside from the stress continuity con- 
ditions (19) no conditions of symmetry were assumed with regard to the 
crack geometry and the external loads a(x-j) and t(x-|). Thus, the inte' 
gral equations (26) derived in Section 2 are valid basically for any 
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number of coll inear cracks defined by X 2 = 0, a^ < < bj, (j=l,...,n) 

along the x-j axis with the additional single-valuedness conditions of 
the form (23) for each crack, namely 
b. 

I ® ’ (k=1.2; j=l,...,n) . (35) 

The only change in the integral equations would be in replacing the 
interval (-a, a) by the sum of the intervals = (a^.,b.) (i=l,...,n) 
corresponding to the col linear cracks. 

As an example, we consider the case of two synmetrically located 
and symmetrically loaded col linear cracks. That is, we assume that 

ai“a, b-j=b, a2~~b, b2~“a, Hi-H2~H, 0i2(Xi,O) — T(x-j) - t(~x^), O22(x-j,0) 

o(x-|) = a(-x-j), (a < x-| < b). In this case, referring to (26) it 
may be shown that for the loading conditions a(x-j) j* 0, t(x^) = 0 and 
o(x-j) = 0, t(x.|) 0 the integral equations uncouple and, by using 

the symmetry conditions, may be expressed as 
b 

^ (36) 


U 

\ *^t.j-x^ * t^^" t(x^), a<x^<b. 


(37) 


where 


^ll^^l’^l^ “ 1 ?i' 1 ) “ 1^11 (x-j j”!--] ) » 


^22^^!*^!^ ^22^^1*^1^ ” 9 ^ (x-jit-j) < b . (38) 

The integral equations (36) and (37) are again solved under the following 
single-valuedness conditions: 
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U u 

G,(t,)dt, “ 0 • I = 0 . 


( 39 ) 


In this problem the unknown functions G-j and G2 are of the form 
F,(x-t) Fo(xn) 


v'(b-x-,)(x.j-aj 


Ab-x^)(x-|-aj 


( 40 ) 


and the stress intensity factors may be obtained from 


Em 


k-jCa) = -y- lim i/2{x^-a) G^(x^) 


x^->a 


k2(a) = Em25 lim >^2(x-|-a) G2(x-|) 


x-j->a 


Em-, 


k-j(b) = - -j- lim i^2(b-x|) G^(x^) 


x-j-^b 


k2(b) = - Em26 lim /2(b-Xi) GgCx-j) 


x^->b 


( 41 ) 


5 . Loading the Strip by a Rigid Wedge 

Consider now the problem described in Figures 1 (b) and (c). It 
is assumed that an orthotropic strip containing a symmetrically located 
crack is loaded by a rigid friction less rectangular wedge of thickness 
Vq. Because of synmetry a-|2(x.j,0) = 0 ,(-®<x-j«»), G2(x-j) = 0 , (-a<x-j<a), 
ki2(xi,ti) = 0 = k2i (x-jst^i* (-a<(x-j,t^)<a), the second equation in ( 26 ) 
is identically satisfied, and assuming that the contact is maintained 
along the entire length of the wedge -b<x-j<b, the first integral equation 
may be expressed as 
a 

w I k^^(x^,t^)]G^(t^)dt^ • ^^2) 
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We now note that for the continuous contact problem shown in Fig. 1(b) 
Gl(x-j) = 0, -b<x^<b, a(x-j) = 0, b<|x^I<a , (43) 

and (42) may be expressed as 

-b a 

7 I k^l(x^,t^)]Gl(t^)dt^ + i| k^l(x^,t^)]G^(t^)dt^ 

-a b 

b<|x^I<a . (44) 

By referring to the collinear crack problem described in the previous 
section and by using the conditions of symmetry in loading and geometry 
the integral equation (44) may be simplified as follows: 

a 

^ \ KTl(x^,t^)]G^(t^)dt^ = 0, b<x^<a . (45) 

Equation (45) must be solved under loading condition 
a 

G^(ti)dt^ = -Vq . (46) 

b 

Since the wedge has "sharp" comers in this case, too, the solution 
of (45) has square root singularities at the end points a and b and may 
be expressed as 

Fi (x, ) 

Si(xi) = — ^ - -- . (b<x,<a) (47) 

Aa-x^)(x^-b) ' 

where F-j is a bounded function. 

Intuitively, it is clear that the continuous contact described 
above is valid only for relatively short wedges, as b/a increases separ- 
ation may occur in the center portion of the wedge. This problem may 
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be investigated and the critical ratio of b/a for initiation of separa- 
tion may be determined by calculating the contact pressure along X 2 = 0, 
-b<x-|<b, and by calculating the value of b/a at which this pressure 
becomes zero (presumably at x-j = 0). It must be emphasized that (b/a)^^ 
would be independent of the wedge thickness v^. In (45) and (46) dividing 
the equations by v^ it is seen that S-j/Vq may be determined uniquely. 

Thus, Vq enters into the expression of the contact pressure simply as 
a multiplying factor and would have no effect on its relative distribu- 
tion. Once G-j(Xi) is determined the contact stress may be obtained from 

a 22 (x-,, 0 ) = J [ t^+5^ K^](x^,t^)]G^(t)dt^, 0<x^<b. 

(47) 

For (b/a) > (b/a)^^ the separation would take place along -c<x-|<c, 
X2=0. In this case by defining 

G^Cx^) = gi(x^), b<|xi|<a, G-|(x^) = g 2 (x^), -c<x^<c , (48) 

and again by observing that 

G-j(x^) = 0, |x^l>a, c<lx^|<b , 


a(x-j) = 0, -c<Xi<c, b<lxi|<a , 

the integral equations of the problem may be expressed as 
a c 



* t^+x^ K^l(x^,t^)]gi(t^)dt^ 



(49) 


(Xi jt-j )lg2(t-| )dt-j 


= 0, b<x^<a , 


(50) 



( 51 ) 


c a 

-c b 

K-j 1 (x-j > t-j ) ^Q"! (t-j )dt-] ~ 0 » “C<x-j <c y 

Note that in (50) and (51) only the first integrals have singular ker- 
nels. Also note that at x^ = +c the contact is "smooth". Thus, the 
solution of (50) and (51) is of the following form: 

fi(xi) ■ ^ 

g^(x-|) = _ » 92(^1) “ f2(x-j )/c^-x-|^ (52) 

/(a-x^)(x^-b) 

where f^ and f2 are (unknown) bounded functions having the symmetry pro- 
perties fj(x-|) = -f^-(-x^), (i=l,2). The loading for the homogeneous 
system (50) and (51) is again provided by 
a 

» 

g^(t^)dt^ "''o • 

b 

In this case, too, the contact stress, indeed the stress component a22(x-j,0) 
may be calculated from (50) or (51) by observing that these expressions 
give ct 22 (xi> 0) outside as well as the inside the cracks. 


6. Results 

In the general problem the system of integral equations (26) is 
solved for loading conditions c(x-j) f 0, x(x-j) = 0 and o(x-j) = 0, t(x-j) f 
0 separately. In doing ^so there is no need for specifying the stiffness 
parameter 6 = explicitly as s enters into the expressions of 

the kernels k.^ through Hfi/a only. TTius, if t{x-j) = 0 one may solve 
for G-j/6 and Gg from (26) and obtain F^(Ta)/6 and p2(+a) for given values 
of H-j/H 2, HS/a, and K. Equations (30) would then give k-j and k2/5, and 
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may be obtained after specifying 5. The important point to make here 
is that decreasing 5 for constant H/a has the same magnifying effect on 
the Mode I stress intensity factor as decreasing the ratio H/a for 
constant 5. Similarly, if c(x^) = 0, tCx-|) / 0 (26) may be solved for 
G-| and 626 for given k, 6 H/a, and and (30) would give kg^and-k-jS. - 

In 'this case, too, decreasing 5 would have the same magnifying effect 
on k 2 as decreasing H/a. The second important observation to be made 
is that in the uncoupled problem, i.e. , for H^=H 2 =H the stress intensity 
factors k^ and k 2 depend on 6 through H5/a only. Therefore, in most of 
the results given in this paper, the stiffness ratio 6 is not explicitly 
specified. 

The results for a single crack are given in Figures 2-6 and Tables 
1-4. Figures 2 and 3 show the effect of the crack location H-j/H on 
the stress intensity factors in a uniformly loaded specific orthotropic 
material, namely yellow birch, and in an isotropic strip. As expected, 
for (H-j/H) 0, k-| and k 2 become unbounded, A somewhat unexpected result 

here is that in the case of shear loading as H-j/H decreases k 2 does not 
monotonously increase, it rather goes through a minimum before becoming 
unbounded (Fig. 3). Figures 4 and 5 show the effect of H/a on the stress 
intensity factors in an isotropic strip (5 =k= 1) and in orthotropic strips 
having the properties k= 1, fi=l/3 or k= 1, 6=3. Note that this covers a 
range of l/81<Eii/E22<81 for the stiffness ratio. For constant shear 
parameter < the effect of H 6 /a on the stress intensity factors is shown 
in Fig. 6 , 

Further results for an orthotropic strip containing a single crack 
are given in Tables 1-4. Table 1 shows the stress intensity factors for 
a concentrated wedge loading (see, for example, the insert in Fig. 10). 
The effect of H/a on the stress intensity factors in an isotropic and 
in two specific orthotropic strips containing a symmetric crack under 
various loading conditions is shown in Table 2. The effect of shear 
parameter k is shown in Table 3 and that of 6 (for a given thickness 
ratio H/a = 0.35) is shown in Table 4. It is seen that compared to 
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the influence of 6 the effect of < on the stress intensity factors is 
negligible. 

A sample result for two pressurized coll inear cracks in an ortho- 
tropic strip is shown in Fig, 7 which also shows the stress intensity 
factors in the corresponding infinite plane (i.e., for H/(b-a)/2 = ~) 
obtained from (see, for example, [11]) 

k,(b) = fT47F>’'^ 1^’ - > 

k,U) = p/rH3577?(|-f^)^^Cf[|f-(f)V. (54) 

where k = /1-aVb^ and K(k) and E(k) are the complete elliptic integrals 
of the first and second kind, respectively. Tables 5-7 show further 
results. In these examples the half crack length (b-a)/2 is used as 
the normalizing length parameter, a/b = 0 and a/b = 1 correspond to 
two limiting cases of a single crack of length 2b and b-a, respectively. 
Perhaps the most interesting result here is the fact that for small 

values of H6/(b-a)/2 k-j (particularly k-j(a)) goes through a minimum as 
a/b goes from zero to 1. This reduction is apparently due to the inter- 
action of the stress fields developing around the two cracks. For 
example, as seen from Figures 8-10 the cleavage stress a22(x-|,0) in the 
plane of the crack, even though tensile and unbounded in the close neigh- 
borhood of the crack tip, becomes compressive in a certain interval away 
from the crack tip. This is essentially due to the "bending" effect 
of the two halves of the strip. This effect diminishes and disappears 
as the thickness ratio H/a increases. No such behavior is observed for 
strips under shear loading. 

Fig. 7 and Tables 5-7 show that for a-»0 as expected k^-(a), (i=l,2), 
becomes unbounded, whereas k^.(b) tends to the values corresponding to 
a single crack of length 2b. It may be observed that the rate of 
increase of k^. (b) (i=l,2) as (a/b) approaches zero is very high. In 
fact, the derivative of k.(b) at a/b = 0 may be shown to be infinite. 
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For example, noting that for small values of a/b 
K(k) = logCi^) + I DogC^)-l]Ca/b)^ + ... , 
ECk) 1 + ^ “ ^](a/b)^ .... 


from C54) it can be shown that 


3^k,(b) : - pv^ ’ 


an 2 


F tVog §] 


C55) 


C56) 


which becomes unbounded at a/b = 0. 

The results for the orthotropic strip loaded by a rectangular rigid 
wedge of thickness are shown in Figures 11-14 and Tables 8 and 9. For 
<=2 and H6/a = 0.6 the contact stress distribution a22(x-j,0) = -p(x), 
(x=x-j/a) corresponding to b/a ratios greater than, equal to, and smaller 
than (b/a)^^ is given in Fig. 11, where (b/a)^^ corresponds to the 
initiation of interface separation. It is seen that for the case of 
continuous contact (e.g., b/a = 0.2) the pressure is minimum at x^=0 
where it becomes "negative" for (b/a) > (b/a) . In the case of b > b 

w ■ w I 

the contact stress is of the following form 


/ X *c 


(57) 


where h(x^) is a bounded function. The calculated length 2c of the 
separation zone is shown in Fig. 12. Critical wedge length b^^ corre- 
sponding to k= 2 and various values of Hs/a is given in Table 8. The 
table also shows the normalized crack tip stress intensity factor 
k-|(a) = k-j^ and the stress intensity factor for the contact stress 
at Xi=b defined by (see (57)) 

kiB = k^Cb) = lim /F-x^ O22(x^,0) . (58) 
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Table 9 shows and k^g for the case of continuous contact, i.e., for 
b<bcr» Additional results for the separation case are given in Figures 
13 and 14. In these figures the starting points of the graphs corre- 
spond to Ch/a)j,^. 


References 

1. S. Krenk, Int. J. Solids Structures, Vol. 11 0975), 449. 

2. D.D. Ang and M.L. Williams, Journal of Applied Mechanics, Vol, 28, 
Trans. ASME 0961), 372. 

3. G.C. Sih, P.C. Paris, G.R. Irwin, International Journal of Fracture 
Mechanics, Vol, 1 (1965), 189, 

4. F. Delale and F.Erodgan, Journal of Applied Mechanics, Vol. 44, 
Trans. ASME (1977), 237. 

5. A.C. Kaya and F. Erdogan, International Journal of Fracture, 

Vol. 16 (1980), 171. 

6. F. Delale, I. Bakirtas, and F. Erdogan, Journal of Applied 
Mechanics, Vol. 46, Trans. ASME (1979), 90.- 

7. S.G. Lekhnitskii, Theory of Elasticity of an Anisotropic Elastic 
Body, Holden-Day Inc. (1963). 

8. S. Krenk, Journal of Composite Materials, Vol. 13 (1979), 108. 

9. A. Cinar, The Crack and Crack-Contact Problems for an Orthotropic 
Strip , M.S. Thesis, Lehigh University, (1980). 

10. F. Erdogan, Mechanics Today , Vol. 4, S. Nemat-Nasser, ed., Pergamon, 
Oxford (1978TrT: 

11. F. Erdogan, Proc. 4th U.S. National Congress of Applied Mechanics, 
ASME (1962), 547, 


- 18 - 



Table 1. Stress intensity factors in an orthotropic strip subjected 
to concentrated forces on the crack surfaces, a 22 (xi» 0 ) = 

Ps(x^) a.j2(x^»0) = 0, K = 1.2895. 



H^=H2=H 

H^=0,7H 

H^=0.4H 

6H/a 

h 

•^1 

-k2 

•^1 

-k2 

(P/a)^ 

(P/a)i/a 

(P/a)y^ 

(P/a)»^ 

(P/a)»^ 

1.5 . 

0.613 

0.681 

0.104 

1.061 

0.353 

1.0 

0.880 

1.008 

0.207 

1.638 

0.700 

0.6 

1.517 

1.767 

0.472 

3.038 

1.555 

0.4 

2.470 

2.895 

0.869 

5.164 

2.875 

0.25 

4.514 

5.344 

1.777 

9.872 

5.864 

0.15 

9.091 

10.822 

3.859 

17.866 

11.193 

0.10 

16.177 

19.098 

7.069 







Table 2. Stress intensity factors in isotropic and orthotropic 
strips which contain a symnetrically located crack 
(i.e,, H]=H2=H) subjected to various loading conditions. 


H/a 





=.q 

k^/p*^ 

k^/{P/a)/a 

kg/qv^ 

<o 

II 

It 

«o 

II 

OJ 

II 

k =2, 6=1/2 

k =1=5 

“O 

II 

CM 

II 

a 

ic =2. 6=1/2 

k =1=6 

k =2=6 

k =2, 6=1/2 

0.2 

9.376 

4.047 


6.249 

2.403 


2.629 

1.789 

. 

0.4 

4.146 

2.175 

8.865 

2.511 

1.133 

5.811 

1.934 

1.391 

2.405 

0.6 

2.758 

1.655 

5.509 

1.536 

0.774 

3.421 

1.640 

1.238 

2.015 

0.8 

2.149 

1.422 

4.047 

1.128 

0.612 

2.403 

1.473 

1.159 

1.789 

1.0 

1.816 

1.294 

3.248 

0.892 

0.524 

1.849 

1.365 

1.113 

1.640 

1.2 

1.611 

1.216 

2.751 

0.745 

0.470 

1.508 

1.289 

1.084 

1.533 

1.6 

1.377 

1.130, 

. 

2.175 

0.581 

0.411 

1.133 

1.192 

1.051 

1.391 

3.0 

1.122 

1.040 

1.466 

0.404 

0.351 

0.643 

1.067 

1.016 

1.175 



Table. 3: The effect of the shear parameter k on the stress intensity factors, 
(a) ^=1.0, ^= 0.35 (b) = 0.4, ^ = 0.35 



a(Xi)=P 

t(x^)=0 

aCx^)=P6(x^) 

t{x^)=0 

a(x-,)=0 

T(x^)=q 

K 


’^2A""'^2B 

CD 

II 

< 

r— 

*^2A‘""'^2B 

*^1A^“'^1B 

k2A=l<2B 

he 

*^2B 

hi 

•^2B 

fiklB 

^2B ' 


pi^ 

6p/a 

(P/a)v^ 

6(P/a)»^ 

q /a 

q/a 

1 

9.519 

-4.855 

6.342 

-3.646 


IS 

2 

8.986 

-4.404 

5.903 

-3.295 


mm 

4- 

8.447 

-3.908 

5.432 

-2.905 


1.441 

8 

7.968 

-3.370 

4.977 

-2.487 

0.563 

1.336 

12 

7.753 

-3.095 

4.737 

-2.266 

0.508 

1.283 



o(x^)=p 

t(x^)=0 

a(x^)=P6(X|) 

t(x^)=0 

a(Xi)=0 

r(x^)=q 


hA^'hB 

hA'^^lB 

l<2A=k2B 

1C 

he 

ho 

*^2B 


p/a 

(P/a)/a 

q/a 

1 

B 

2.971 

2.047 

2 


2.826 

1 .890 

4 


2.697 

1.717 

8 

wm 

2.612 

1.550 

12 

4.611 

2.592 

1.463 

16 

4.692 

2.613 

1.405 











Table 4. 


Effect of the -stiffness parameter 5 
on the stress intensity factors. 

Hi H2 

— = -:r = 0.35. < = 2. 


6 


^22^^1 ’0)~P5(x-j ) 


^li 

p/a 

(P/a)*^ 

'^2B 

q/a 

0.3 

20.350 

14.218 

3.231 

0.4 

13.880 

9.426 

, 2.829 

0.6 

8.348 

5.436 

2.353 

0.8 

5.954 

3.732 

2.075 

1.2 

3.851 

2.261 • 

1.755 

1.5 

3.103 

1.745 

1.610 

1.8 

2.637 

1.427 

1.507 

2 

2.415 

1.279 

1.453 

3 

1.797 

0.858 

1.282 

10 

1.112 

0.393 

1.043 
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Table 5. Stress intensity factors in an orthotropic 

strip containing two identical collinear cracks 
loaded by unifom crack surface pressure p 
or shear q; Hi=H 2 =H, k= 1, H6/(b-a)/2 = 0.4. 



<^22^^1 

.o)=-p 

ai2(Xi 

.0)=-q ■“ 

2a 

k^(b) 

k^(a) 

l<2{b) 

kgCa) 

b-a 





0 

9.376 

00 

2.629 

00 

.01 

3.693 

6.996 

2.106 

5.837 

.1 

3.788 

2.837 

1.952 

2.300 


3.962 

3.113 

1.935 

1.989 


4.074 

3.642 

1.933 

1.939 


4.124 

3.971 

1 .933 

1.933 


4.138 

4.103 

1.933 

1.932 


4.141 

4.138 

1.933 

1.933 


4.140 

4.143 

1.933 

1.933 

.8 

4.140 

4.142 

1.933 

1.933 

.9 

4.139 

4.140 

1.933 

1.933 

1 

4.139 

4.140 

1.933 

1.933 


4.142 

4.142 

1.933 

1.933 
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Table 6. Stress intensity factors in an orthotropic 

strip containing two identical collinear cracks 
loaded by uniform crack surface pressure p 
or shear q; <=2, H6/(b-a)/2 = 0.4. 


2a 

c^22 ^ 

Q. 

II 

O 

^22 ^^1 

0 

II 

1 










q(^) 

0 

8.865 

CO 

2.405 

00 

.01 

3.626 

7.023 

1.993 

5.657 

.1 

3.697 

2.890 

1.825 

2.252 

.2 

3.847 

3.126 

1.799 

1.908 

.3 

3.950 

3.547 

1.792 

1.824 

.4 

4.004 

3.814 

1.790 

1.800 

.5 

4.028 

3.947 

1.790 

1.793 

.6 

4.039 

4.005 

1.790 

1.790 

.7 

’4.044 

4.030 

1.789 

1.790 

.8 

4.045 

4.040 

1.789 

1.790 

.9 

4.046 

4.044 

1.789 

1.789 

1 

4.046 

4.046 

1.789 

1.789 

2 

4.047 

4.047 

1.789 

1.789 
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Table 7. Stress intensity factors in an orthotropic 

strip containing two identical col linear cracks 
loaded by uniform crack surface pressure p 
or shear q; k= 2, H(5/(b-a)/2 = 0.8. 




Q. 

1 

II 

O 

Q 

ro 

X 

o 

II 

2a 

k,(b) 

k^(a) 

kgCb) 

kgCa) 

b-a 

P(^)^ 



q(^) 

0 

4.047 

00 

1.789 

00 

.01 

2.207 

6.843 ■ 

1.677 

4.951 

.1 

2.012 

2.448 

1.482 

2.057 

.2 

2.017 

1.999 

1.433 

1.676 

.3 

2.046 

1.913 

1.412 

1.532 

.4 

2.079 

1.931 

1.402 

1.464 

,5 

2.106 

1.981 

1.397 

1.430 

.6 

2.127 

2.033 

1.394 

1.412 

.7 

2.143 

2.075 

1.393 

1.402 

.8 

2.154 

2.107 

1.392 

1.397 

.9 

2.161 

2.130 

1.391 

1.394 

1 

2.166 

. 2.146 

1.391 

1.393 


2.174 

2.174 

1.391 

1.391 


2.175 

2.175 

1.391 

1.391 
















Table s. The critical values of ^ for 

a 

various values of 6H/a, shear para- 
meter K = 2. 


6H 

a 

(|> 

cr 

6a “^lA 

6a *^18 

1.5 

0.961 

4.761x10“^ 

-4.568x10"^ 

1.0 

.0.665 

1.569x10"^ 

-1.476x10"^ 

0.6a 

0.293 

8.544x10"^ 

-6.7,92x10"^ 

0.40 

0.161 

5.836x10"^. 

-3.874x10“^ 

0.20 

0.059 

2.801x10"^ 

-1.188x10'^ 
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Table 9. Stress intensity factors in an orthotropic 
strip loaded by a rectangular wedge for 
(b/ai< (b/a)cr (Fig. lb), k=2, Icq = 
Evo*WTa5) (underlined values correspond to 
initiation of separation at = 0). 
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Hi/H 

Figure 2: The effect of the crack location on the stress inten- 


sity factors for uniform surface pressure. H * 0.75a, 


6 = 1 = K for the isotropic materials and 6 = 1.1175, 

K =1.2895 for the orthotropic material (yellow birch). 






crack surface. 



Figure 4: Effect of the crack length on the stress intensity 

factor for a symmetrically located crack under uniform 
pressure, k = 1. 






Figure 6: The effect of 6H/a on the stress intensity factors 


for a symmetrically located crack under uniform 
pressure or uniform shear, k = 1.2895. 



5 



Figure 7. Stress intensity factors 
orthotropic strip. 


-o)/2 


)/p7(b-a)/2 


8H/(b-a)/2=0.4 



■ay/2 

H = cd 

■a)/2 


I .1 

0.3 0.4 


for two coll inear cracks in an 





I 2 3 4 5 6 


x,/a 

Figure 8. The effect of the relative crack size on the normal stress 
^22 the plane of the crack in an orthotropic strip con- 
taining a pressurized crack (see the insert in Fig. 9 for 
notation). 



Figure 9. The effect of material orthotropy on the normal stress 
a 22 (x-|, 0 ) in a strip containing a pressurized crack (6 
K=1 isotropic strip). 




1.0 



Figure 10. Same as Figure 8 in a cracked strip under concentrated 
wedge loading. 




Figure 11. Pressure distribution for the wedge problem, 5H = 0.6a, 
k=2, p(x) = -022^^1’^^’ x=x^/a. 






Figure 12: Separation length for the wedge problem, k = 2. 




Figure 13: Crack tip stress intensity factor for the wedge prob- 


lem, K = 2. 




Figure 14: Wedge tip stress intensity factor for the crack con- 
tact problem, k = 2. 
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